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Abstract Wc consider an optimal control problem describing a laser-induced population transfer 
on a n-level quantum system. 

For a convex cost depending only on the moduli of controls (i.e. the lasers intensities), we prove 
that there always exists a minimizer in resonance. This permits to justify some strategies used 
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in experimental physics. It is also quite important because it permits to reduce remarkably the 
complexity of the problem (and extend some of our previous results for n — 2 and n = 3): instead 
£C) , of looking for minimizers on the sphere 5* 2 "" 1 C C" one is reduced to look just for minimizers on 

the sphere S" 1 " 1 C K". 

Moreover, for the reduced problem, we investigate on the question of existence of strict abnormal 
minimizer. 
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1 Introduction 

The problem of designing an efficient transfer of population between different atomic or molecular levels is 
crucial in many atomic-physics projects [17, 26, 20, 35, 6]. Often excitation or ionization is accomplished by 
using a sequence of laser pulses to drive transitions from each state to the next state. The transfer should be as 
efficient as possible in order to minimize the effects of relaxation or decoherence that are always present. In the 
recent past years, people started to approach the design of laser pulses by using Geometric Control Techniques 
(see for instance [3, 16, 32, 24, 25, 33]). Finite dimensional closed quantum systems are in fact left invariant 
control systems on SU(n), or on the corresponding Hilbcrt sphere S* 2 " -1 c C", where n is the number of 
atomic or molecular levels. For these kinds of systems very powerful techniques were developed both for what 
concerns controllability [4, 18, 22, 23, 34] and optimal control [1, 21]. 

The most important and powerful tool for the study of optimal trajectories is the well known Pontryagin 
Maximum Principle (in the following PMP, see for instance [1, 21, 31]). It is a first order necessary condition for 
optimality and generalizes the Weierstrafi conditions of Calculus of Variations to problems with non-holonomic 
constraints. For each optimal trajectory, the PMP provides a lift to the cotangent bundle that is a solution to 
a suitable pseudo-Hamiltonian system. 

Anyway, giving a complete solution to an optimization problem (that for us means to give an optimal 
synthesis, see for instance [7, 9, 14, 15, 30]) remains extremely difficult for several reasons. First one is 
faced with the problem of integrating a Hamiltonian system (that gencrically is not integrablc excepted for 
very special costs). Second one should manage with "non Hamiltonian solutions" of the PMP, the so called 
abnormal extremals . Finally, even if one is able to find all the solutions of the PMP it remains the problem 
of selecting among them the optimal trajectories. For these reasons, usually, one can hope to find a complete 
solution of an optimal control problem in low dimension only. 

1 This research has been supported by a Marie Curie Fellowship. Program "Improving Human Research Potential and the 
Socio-economic Knowledge Base", Contract Number: HPMF-CT-2001-01479. 

Supported by european networks NCN (ERBFMRXCT97-0137) and NACQ2 (N.HPRN-CT-1999-00046) 
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This paper is the continuation of a series of papers on optimal control of finite dimensional quantum systems 
[10, 11]. The main purpose is to show that for a certain class of quantum systems (that contains systems useful 
for applications) one can reduce remarkably the complexity of the problem. More precisely we prove that for 
a convex cost depending only on the moduli of controls (e.g. amplitude of the lasers): 

• there always exists a minimizer in resonance that connects a source and a target defined by conditions 
on the moduli of the components of the wave function (e.g. two eigenstates, see Theorem 1). Roughly 
speaking to be in resonance means to use lasers oscillating with a frequency equal to the difference 
of energy between the levels that the laser is coupling. As a consequence one gets a reduction of the 
dimension of the problem from 2n — 1 to n — 1 (n being the number of energy levels), see Corollary 1. 
From a physical point of view this means that one is reduced to look just for the amplitudes of the lasers; 

• for the reduced system in dimension n — 1, we prove that close to any time t of the domain of a given 
minimizer, there exists an interval of time where the minimizer is not strictly abnormal (see Theorem 
3). This result is a first step in trying to prove the conjecture that for the reduced problem there are no 
strictly abnormal minimizer. 

This extends some of our previous results (see [10, 11]). 

Here we are considering a class of systems on which it is possible to eliminate the so called drift term. 
This includes n-level quantum systems in the rotating wave function approximation (RWA) and in which each 
laser couples only close levels. For this kind of systems our reduction is crucial to give a complete solution 
for systems with n = 3 for several costs interesting for applications. Moreover it gives some hope to find the 
time optimal synthesis for systems with four levels and bounded lasers. Finally it is of great help in finding 
numerical solutions to problems with n > 4. 

The paper is organized as follows. In Section 1.1 we introduce the physical model. In order to characterize 
controllability, in Section 1.2 we associate a topological graph to the system in a unique way. In Section 1.3 
we formulate the optimal control problem, while in Section 1.4 we discuss the costs that are interesting for 
applications. In Section 1.5 we introduce some key definitions and we state our main questions. 

Section 2 is devoted to recall how to eliminate the drift term from the control system, by means of a unitary 
time-dependent change of coordinates plus a unitary change of controls (interaction picture) . This permits to 
write the system in "distributional" form. 

Our main results are stated in Section 3 (see Theorems 1,3). Beside stating that there are always minimizers 
in resonance and studying strictly abnormal minimizers for the reduced problem, we investigate also the question 
if every minimizer is in resonance (see Theorem 2). More precisely, we state that, under the assumption that the 
cost function is strictly increasing with respect to the moduli of the controls, every minimizer is weakly-resonant 
in a suitable sense. 

In Section 4 we prove the main results about resonance. To investigate resonance, the key point is to 
identify the components of the controls responsible of the evolution of the moduli of the components of the 
wave function. The difficulty comes from the fact that, when a coordinate is zero, the dimension of the 
"distribution" may fall. This obliges to divide the domain of a minimizer in suitable sub-domains (see Section 
1.5.1). In the first part of Section 4 we prove our main results, while in Section 4.1 we give an alternative 
geometrical interpretation. In Section 4.2 we investigate the question if it is possible to join every couple of 
points in S" 2 "" 1 C C™ by a resonant minimizer. In Section 4.3 we give an example of a minimizer (for a non 
strictly increasing cost) which is not weakly-resonant and we propose an open question. 

The most technical part is Section 5 where we investigate about strictly abnormal minimizers for the reduced 
problem on the sphere S"^ 1 C K". For trajectories whose coordinates are never zero, our results are just a 
consequence of the fact that the dimension of the distribution is the same as the dimension of the state space. 
But again the difficulty is coming from the fall of the dimension when some coordinates are zero. 

1.1 The Model 

In this paper we consider closed finite dimensional quantum systems in the rotating wave function approximation 
(RWA). More precisely we consider systems whose dynamics is governed by the time dependent Schrodingcr 
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equation (in a system of units such that h = 1): 

at ^ 

4>(.) := (M-)> ■■■>»(■)) : R — » C", £. l^l 2 = 1, 
where it holds: 

(HI) D = diag(Ei, ...,E n ) and V(t) is an Hermitian matrix (V(i)j t fc = V'(t)^-), measurable as function of t, 
whose elements V{t)j^ are either identically zero or controls . 

Here (*) indicates the complex conjugation involution. E\, E n are the energy levels of the quantum system 
and the controls V^fc(.), that we assume to be C-valucd measurable functions, are different from zero only in a 
fixed interval [0, T]. They are connected to the physical parameters by Vj.k(t) — Hj t kJ 7 j i k(t)/2, j, k = 1, ...,n, 
with Tj t k the external pulsed fields (the lasers) and fij.k = fik.j > the couplings (intrinsic to the quantum 
system). In the following we say that two levels E^E^ are coupled if Vj^ is a control (and not zero). Moreover 
if all the couplings fj,j t k are equal to a constant (i (that we normalize to 1) we say that the system is isotropic. 
Otherwise we say that the system is non-isotropic. 

Remark 1 The term D = diag(Ei, ...,E n ) in equation (1) is called drift and it will be eliminated in Section 
2, thanks to the fact that we are in the rotating wave approximation and that each control couples only two 
levels. 

Remark 2 This finite-dimensional problem can be thought as the reduction of an infinite-dimensional problem 
in the following way. We start with a Hamiltonian which is the sum of a "drift-term" D, plus a time dependent 
potential V(t) (the control term, i.e. the lasers). The drift term is assumed to be diagonal, with eigenvalues 
(energy levels) ... > Es > E2 > E±. Then in this spectral resolution of D, we assume the control term V(t) to 
couple only a finite number of energy levels by pairs. Let {Ej}j e j (I finite) be the set of coupled levels. The 
projected problem in the cigcnspaccs corresponding to {Ej}j^i is completely decoupled and it is described by 
(1), (HI). 

The function ip : [0,T] — > C C" solution of the Schrodinger equation (1) is called the wave function . 

The physical meaning of its components tpjif) is the following. For time t < and t > T (where the controls 
are zero), \i/jj(t)\ 2 is the probability of measuring energy Ej. Notice that, in the intervals of time where V is 
identically zero, we have: 

^■(t)| 2 = 0, j = l,...n. 
In the following a state for which we have = 1 for some j will be called an eigenstate. 

Remark 3 This model is physically reasonable in the case in which: 

• the number of energy levels is not too big, and they are distinct by pairs; 

• there are not too many couplings between the energy levels. 

The most interesting case is perhaps the one in which only close levels are coupled. In this case V has non null 
coefficients Vj t k only on the second diagonals (for j and k such that \j — k\ = 1) and the Hamiltonian reads: 

/ Ei V ia {t) ••• \ 

v^t) e 2 v 2 , 3 (t) ■■. ; 

V 2 * 3 (t) ■■■ ■■■ (2) 



H(t) 



'■ '■ '■ E n -\ V n -i : n(t) 

V ••• V*_ x Jt) E n 



In the rest of the paper, we will refer to this model as to the most important example. 
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1.2 Controllability and Representation of the System with a Graph 



To each system of the form (1), (HI) one can associate a topological graph (i.e. a set of points and a set of edges 
connecting the points) in a very natural way. The points are associated to the energy levels Ej (j = 1, ...,n) 
and two points Ej, Ek are connected by an edge iff the clement Vj.k is a control. This makes sense also if all 
the energy levels are zero (it happens after elimination of the drift, see Section 2). 

Lifting the problem for the operator of temporal evolution (i.e. on U(n)) and with standard arguments of 
controllability on compact Lie groups and corresponding homogeneous spaces, one gets the following (see the 
recent survey [34] or the papers [4, 10, 11, 13, 18, 22, 23]): 

Proposition 1 The control system (1), under the assumption (HI) is completely controllable from any initial 
to any final condition if and only if the corresponding graph is connected. 

In the following we deal with optimal control problems and we assume existence of minimizers for every couple 
of points. Hence we assume: 

(H2) The graph associated with the control system (1),(H1) is connected. 

If (H2) does not hold, then all the results of the paper are true for the restricted systems associated to the 
connected parts of the Graph. 

Remark 4 Notice that to guarantee controllability it is necessary (but not sufficient) to have at least n — 1 
controls Vj.k (with j < k). 

1.3 The Optimal Control Problem 

In this paper, we are faced with the problem of finding optimal trajectories for a convex cost depending only on 
the moduli of controls between a source and a target defined by conditions on the moduli of the components 
of the wave function. These are the most common situations in physics since the squares of moduli of the 
components of the wave function represent the probabilities of measures. More precisely our problem is the 
following: 

Problem (P) Consider the control system (1),(H1),(H2) and assume that for time t = the state of the 
system is described by a wave function ip(0) whose components satisfy (|^i(0)| 2 , \ip n (0)\ 2 ) G Si n , where Si n 
is a subset of the set: 

S := {(oi, ....a*) G (R+y i : a * = X >- ( 3 ) 

We want to determine suitable controls Vj.k(-), j,k = l,...,n, defined on an interval [0,T], such that for time 
t = T, the system is described by the wave function ip(T) satisfying (\ipi(T)\ 2 , |-0„(T)| 2 ) G Sfi n , where 
Sfi n C <£>, requiring that these controls minimize a cost that depends only on the moduli of controls: 



I f(V(t))dt. (4) 
Jo 

Remark 5 Typical sources and targets are eigenstates. For instance if the source and the target are respectively 
the eigenstates 1 and n, we have Si n = (1, 0, 0), Sfi„ = (0, 0, 1). 

In the following to guarantee the existence of minimizers we assume: 

(H3) The function f° (that depends only on the moduli of controls) is convex. Si n and Sfi n are closed 
subset of iS>. Moreover the moduli of controls are bounded: there exist constants Mj k > such that 
\V jtk (t)\ < M jtk , for every t G [0,T]. 

Remark 6 In the problem (P), the final time can be fixed or free depending on the problem and on the explicit 
form of f° (see below). 
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Remark 7 Notice that if f° is convex and depends only on the moduli of controls then it is an increasing 
function (as function of the moduli of controls). Hypotheses (H3) could be relaxed, in particular if we do not 
require existence of minimizcrs for each couple of points, it is not always necessary to assume that controls 
are bounded or that /° is convex. Moreover the hypotheses of boundedness of controls could be changed with 
suitable hypotheses of growing of / at infinity. Anyway these investigations are not the purpose of the paper. 
The costs on which we are interested are those described in the next Section, that are convex (some of them 
strictly). Moreover the corresponding minimization problems are always equivalent to minimization problems 
with bounded controls. 

In the following if V(.) is a function satisfying (HI), (H2), (H3), and ip(.) the corresponding absolutely contin- 
uous trajectory we say that the couple V(.)) is an admissible pair. 



1.4 The interesting costs 

In this paper, we consider only costs that depend on the moduli of controls, since these are the interesting costs 
from a physical point of view, like those described in the following. 



Energy 



fdt, f = j2-^\v jtk \ 



(5) 



This cost is proportional to the energy of the laser pulses. After elimination of the drift (see Section 2), the 
problem (P), with this cost is a sub-Riemannian problem or a singular-Ricmannian problem depending on the 
number of controls (see (HI)) and was studied in [10, 11], in the isotropic case (fJLj } k = n), for n — 2,3 and 
in the case in which the Hamiltonian is given by (2). For sub-Riemannian and singular-Ricmannian geometry 
see for instance [5, 19, 27]. y For this cost the final time must be fixed otherwise there are not solutions to 
the minimization problem. We recall that a minimizer for this cost is parametrized with constant velocity 



<k 



V 



j,k\ 



const) and it is also a minimizer for the sub-Riemannian length: 



fdt, f 



(0) 



This cost (6) is invariant by re-parametrization, hence the final time T can be equivalcntly fixed or free. For 
the costs (5), (6) the controls arc assumed to be unbounded . Anyway if the final time T is fixed in such a way 



the minimizer is parametrized by arc-length (^2 



j<k 



\Vj,k\ — 1); then minimizing the cost (5) is equivalent 



to minimize time with moduli of controls constrained in the closed set: 



^ u 2 



V 



j,k\ 



< 1. 



(7) 



In the isotropic case 
isotropic-length. They read: 



fi = 1 the energy and the length are called respectively isotropic-energy and 



[ T fdt, f = 52\V Jt 



If,, 



(8) 
(9) 



and the set (7) becomes a closed ball. In the non-isotropic case we call (5) and (6) respectively non-isotropic-energy 
and non-isotropic-lcngth. 
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Isotropic and Non-isotropic Area 

/ f° dt, f° = \Vj k\, (isotropic case Hi = /i = 1) (10) 

Jo W' fe 

rT 

These costs are proportional to the area of the laser pulses J Q X^<fc l-^i,fc(^)l dt. They are invariant by 
re-parametrization. Minimizing these costs is equivalent to minimize time with the following constraint on 
controls: 

Time with Bounded Controls 

If we want to minimize time having a bound on the maximal amplitudes of the lasers, \!Fj,k\ < v j,k (j < k) 
then we get the bounds on Vj^' \Vj,k\ < A*.j,fc i/ i,fe/2- With the change of notation fij.k^j^/^ — > Hj,k we get: 

\Vj,k\<H>k- (11) 

In this case we have f° = 1. Now the isotropic case [ij,k = fi = 1 corresponds to an isotropic system in which 
all the lasers have the same bound on the amplitude. 

In this case the problem is equivalent (up to a re-parametrization) to minimize: 



/ f° dt, f° = max 



t-T 

(12) 



In the following figure the level sets /° = 1 for the n — 3 case for the Hamiltonian (2) are drawn (of course in 
M 2 instead of C 2 ). If we consider the equivalent problems of minimizing time with constrained controls, then 
controls should lie in the compact regions whose borders are these level sets. 
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Isotropic Energy: f° = | Vx,2 1 2 + 1^2,3 1 2 



Non-isotropic Energy: f° = + 

Ml, 2 £^2,3 








V 


2,3 











Isotropic Area: / = |T4,2| + 1^2,3 



Non-isotropic Area: / 



_ Vl,2| _|_ |V2,3| 
Ml, 2 A*2,3 










' ^2,3 








/ ^1,2 


max(|Vi )2 | 


.[^2,3|) 


Time with b. 


contr: f° = 



■\ViM |v 2 , 3 | 



Ml, 2 ' M2,3 



(Isotropic Case) 





1 V2,3 


















v 1>2 









(Non-isotropic Case) 

^2,3 



Vi:. 



Figure 1: level /° = 1 for different cost functions 
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1.5 Main Questions 



Although the system (1), (HI) has a lot of symmetries, and good properties, the problem (P) in general is 
very difficult to solve also for the costs described above, for the reasons explained in the introduction: 

• problem of integrability of the Hamiltonian system associated with the PMP, 

• presence of abnormal extremals, 

• problem of selecting the optimal trajectories among those satisfying the PMP. 

Up to now the problem is solved for n — 2 (for the costs described in Section 1.4, that, in this case, are all 
equivalent) and for n = 3, for the Hamiltonian given by (2) for the isotropic energy (see [10, 11]). In both 
cases optimal controls appear to be in resonance (with the difference between the energy levels that they are 
coupling) and every strictly abnormal extremal is not optimal. In this paper we generalize these results for 
more general systems and costs. Let us first introduce the definition of resonance and abnormal extremals: 

Definition 1 (Resonance) Consider the control system (1),(H1) and an admissible pair (ip(-), V(.)) defined 
in an interval [0,T]. We say that the couple (ip(.),V(.)) is in resonance (or is resonant ) if the controls Vj t k(-) 
have a.e. the form: 



In formula (15), ifipj(0) = 0, then arg(tpj(0)) is intended to be an arbitrary number. 

Remark 8 From a physical point of view, to be in resonance means to use lasers (described by complex 
functions) oscillating with a frequency given by us/(2ir) where us is the difference of energy between the levels 
that the laser is coupling. Notice that in formula (14) Uj^i-) takes real values (it describes the amplitude of 
the lasers). 

In the definition of resonance, the phases (pj t k sometimes are important and sometimes not, depending on 
the explicit form of the control term V(t) in (1) and on the choice of the source and the target. For instance 
if we are considering the problem described by the Hamiltonian (2) and we start from an eigenstate (e.g. 
|V»i(0)| = 1) then all the phases tpj^k appearing in formula (13) are arbitrary. In fact in this case we have n — 1 
arbitrary numbers arg(ipj(0)) and n — 1 controls. Therefore in this case it is not necessary to synchronize the 
phases of the lasers between them. 

On the other side, if we start from a state that is not an eigenstate or we have more controls than the mini- 
mum number necessary to guarantee controllability (see Remark 4), to be in resonance we need to synchronize 
the lasers according to formula (15). 

In any case notice that a global factor of phase in front of ip is not important because it can be eliminated 
with a unitary transformation. Hence all the phases tpj.k can be shifted by an arbitrary factor a. 

Remark 9 To prove that a minimizer corresponds to controls in resonance, is very important because, as we 
will see after elimination of the drift, this permits us to reduce the dimension of the problem from the complex 
sphere S' 2 ™- 1 C C™ to the real sphere S*™" 1 C E™ (See the Corollary 1 in Section 3). This reduction of 
dimension is crucial in finding complete solutions to the optimal control problem in many cases. For instance 
for the class of problem described by the Hamiltonian (2): 

A. an isotropic or non-isotropic minimum time problem with bounded controls for a 3Tevel system, is a 
problem in dimension 5. In this case, since the dimension of the state space is big, the problem of finding 
extremals and selecting optimal trajectories can be extremely hard. The possibility of proving that one 
can restrict to minimizers that are in resonance (and this is the case!) permits us to reduce the problem 
to a bi-dimensional problem, that can be solved with standard techniques (see for instance [9, 21, 29, 36]); 



V jtk (t) = Uj^ty^-^^' 2 ^^, where: 

U hk {.) : [0,71 -R, U jtk = -U hJ , 

ip-i.k := arg(i/jj(Q)) - arg{ifs k (0)) G [-it, it] 



(13) 
(14) 
(15) 
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B a nonisotropic minimum energy problem for a 3-level system is naturally lifted to a left invariant sub- 
Riemannian problem on the group SU(3). This problem cannot be solved with the techniques used in 
[11] for the isotropic case, because now the cost is built with a "deformed Killing form". Anyway if one 
can restrict to minimizers that are in resonance, the problem is reduced to a contact sub-Riemannian 
problem on SO(3), that does not have abnormal extremals (since it is contact) and the corresponding 
Hamiltonian system is completely integrable (since it leads to a left invariant Hamiltonian system on a 
Lie group of dimension 3). 

Items A. and B. are studied in a forthcoming paper. 

C. An isotropic or nonisotropic minimum energy problem for a 4-level system is a problem on S 7 or can be lifted 

to a left invariant sub-Riemannian problem on SU (4). In this case we conjecture that the Hamiltonian system 
given by the PMP is not Liouvillc integrable both before and after the reduction to real variables. Any- 
way this reduction can be crucial in looking for numerical solution to the PMP (it is of course easier making 
numerical computation on S 3 than on «S 7 !). See [12] for some numerical solutions to this problem. 

Moreover proving that there always exists a minimizer in resonance permits us to justify some strategies used 
in experimental physics. 

To introduce the concept of abnormal extremals, we have to state the PMP, that is a necessary condition for 
optimality. The proof can be found for instance in [1, 21, 31]. 

Theorem (Pontryagin Maximum Principle) Consider a control system of the form x = f(x,u) with 
a cost of the form f°(x,u) dt, and initial and final conditions given by x{§) G M;„, x(T) G Mfi n , where 
x belongs to a manifold M and u G U C K m . Assume moreover that M, /, f° are smooth and that Mi n 
and Mfi n cire smooth submanifolds of M. If the couple (u(.),x(.)) : [0,T] C M — > U x M is optimal, 
then there exists a never vanishing Held of covectors along x(.), that is an absolutely continuous function 
(P(.),p ) : t £ [0, T] — > (P(t),p ) G T*, t) M x R (where p < is a constant) such that: 

i) x(t)=^(x(t),P(t),u(t)), 

ii) P(t) = -^L(x(t),P(t),u(t)), 
where by definition: 

W(x,P,u) :=< Pf(x,u) > +p f°(x,u). (16) 

Moreover: 

iii) M(x(t),P(t),u(t)) =W. M (x(t),P(t)), for a.e. te [0,T], 

where P) max„ e [/ W(x, P, v). 

iiii) Mj\/(i(t), P(t)) — k > 0, where k depends on the final time (if it is fixed) or k = if it is free. 

v) < P(0),T x (QjMi n >=< P(T),T x (T)Mfi n >= (transversality conditions). 

The real-valued map on T*M x [/, defined in (16) is called PMP-Hamiltonian . A trajectory x(.) (resp. 
a couple (u(.),x(.))) satisfying conditions i), ii), iii) and iiii) is called an extremal (resp. extremal pair.) If 
x(.) satisfies i), ii), iii) and iiii) with po = (resp. po < 0), then it is called an abnormal extremal (resp. a 
normal extremal ). It may happen that an extremal a;(.) is both normal and abnormal. So it makes sense to 
speak of strictly abnormal extremals (extremals that are abnormal but not normal). 

Remark 10 Notice that the definition of abnormal extremal does not depend on the cost but only on the 
dynamics (in fact if po = 0, the cost disappears in (16)). After elimination of the drift (see Section 2) our 
control system (1) will be transformed into a distributional system, that is a system of the form x = ^ UjFj(x). 
For these systems abnormal extremals are also singularities of the end point mapping, and have very special 
features (see for instance the recent monograph [8]). 
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Remark 11 In this paper we are dealing only with non-state dependent costs, i.e. f°(x, u) = f°(u). 

Consider the problem (P) under the assumption (H3). Hence a minimizer exists. The main purpose of this 
paper is to answer to the following questions: 

Ql Does there exist a minimizer that corresponds to controls in resonance? 

Q2 Once we have given a positive answer to question Ql, we would like to answer the question: are all the 
minimizers of the problem (P) in resonance? 

Q3 Do there exist strictly abnormal minimizers? 

The answers to questions Ql is yes, while Q3 is still an open question, although we are able to prove a result 
that strongly restricts the set of candidates strictly abnormal minimizers. 
These results are formalized in Theorems in Section 3. 

The answer to question Q2 is in general no, but for an increasing cost we show that every minimizer is 
in weak-resonance (again this will be formalized in a Theorem after elimination of the drift). To define this 
notion, in the next Section we introduce some notations that will be also useful in Section 4 to answer Ql and 
Q2. 

Finally another interesting question for physicists is: 

Ql' is it possible to join two arbitrary states ip 1 an d ip 2 by an optimal trajectory in resonance? 

The answer is, of course, in general, no. But there arc couples of points for which the answer is yes. In 
particular if we consider eigenstates, it is true (see Section 4.2). 

1.5.1 Weak- Resonance 

To introduce this concept, we need some notations. 

Consider an admissible pair (ip(.), V(.)) and define the following subset of ]0,T[: 

Badj.k = {t e]0.T[, s.t. Vj(t) = or ip k (t) = 0}. (17) 

Since ip(.) is continuous, the sets ]0, T[\Badj^ are open. Then they can be expressed as a countable union of 
open maximal intervals: 

]0,T[\Badj,k=--{Jlj,k,i, j,k = l,...,n, I = 1, ...to, where to 6 {0} UNU {oo} . (18) 
i 

In other words Ij lk j are the maximal open intervals on which tpj(t)ip k (t) 0- 

Remark 12 In general Uilj.k,i is smaller than [0, T] since it may happen that il)j(t)tp k {t) = on some interval 
of positive measure. 

Now we are ready to give the following: 

Definition 2 (weak-resonance) Consider the control system (1),(H1) and an admissible pair (ip(.),V(.)) 
defined in an interval [0,T]. We say that the couple (%/}(.), V(.)) is in weak-resonance (or is weakly-resonant) 
if in each interval Ij.k,i (defined above) the controls Vj tk {.) satisfy a.e.: 

^>*li,.*., (*) = U jAl (t)e^- E ^/ 2+ ^'\ where: (19) 
U jtk!l (.):[0,T]^R, U jikt i = -U kJi i (20) 
<Pj,k,i ■= arg(ipj(a jik! i)) - arg(ip k (a jik! i)) e [-7r,7r] (21) 

where Ij, k ,i ='^\a,j,k,U bj,k,l[- In formula (21), if^j{o-j,k,l) — 0j then arg{{ipj (aj tkt i))) is intended to be an arbitrary 
number. 

Remark 13 Roughly speaking a control Vj lk (.) is weakly-resonant if it is resonant in each interval of time in 
which the states that it is coupling (i.e. ipj and ip k ) are different from zero. 

If the cost is not strictly increasing (for instance minimum time with bounded controls) in general there arc 
minimizer that are not in weak-resonance, see Section 4.3 for an example and another open question. Finally 
notice that a resonant minimizer is also weakly-resonant. 
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2 Elimination of the Drift Term 



The great advantage of the model presented in the previous Section (in which each control is complex and 
couples only two levels) is that we can eliminate the drift term D = diag{E\ 1 E n ) from equation (1) (hence 
getting a system in distributional form x = ■ UjFj(x)), simply by using the interaction picture. This is made 
by a unitary change of coordinates and a unitary change of controls. Since the transformation is unitary, Si n , 
Sfin and the moduli of the components of the wave function are invariant. As a consequence the original and 
the transformed systems describe exactly the same population distribution. 

Assume that tf>(t) satisfies the Schrodinger equation (1). Let U(t) be a unitary time dependent matrix and 
set i/j(t) = U(t)ip' (t). Then ip'(t) satisfies the Schrodinger equation: 



with the new Hamiltonian: 



If we choose: 



W =U- 1 HU-iU~ 1 ^-. (22) 
dt y ' 



U{t) = e~ lU \ (23) 

and we recall that H = D + V(t), we get TC' = e lDt V(t)e~ tDt , that is a Hamiltonian whose elements are either 
zero or can be redefined to be controls. Hence the drift is eliminated. Finally dropping the primes and including 
the i in the new Hamiltonian (ij/ — > ip. H := —iH') the Schrodinger equation reads: 

^=H(t)m- (24) 

Here H is skew-Hcrmitian and its elements are either zero or controls. There is no more drift and the dynamics 
is in "distributional" form, i.e. tpj(t) = J2k=i Hj,k(t)ipk(t)- 

The relation between the old controls V{t)j t k (for equation (1)) and the new controls H(t)j : k for equation 
(24) is the following: 

V{t) hk = (ie- iDt H(t)e iDt ). k = H(t) jtk e *[(^-^)*+V2]. (25) 

Remark 14 Notice that the transformation (23) kills also a not-only-diagonal drift. While to kill a time 
dependent drift we need the transformation (see [10]): 

U(t) = e J *o 



3 Statement of the Main Results 

With the transformation described above (formula (23), and following), the problem (P) becomes the problem 
(P') below and the answers to questions Ql, Q2, Q3 are contained in the next Theorems, that we are going 
to prove in Sections 4, 5. 



Problem (P') Consider the minimization problem: 

#(t) 



dt 
mm 



= H(t)1>(t) (26) 



fin [ f (H) dt. (27) 
Jo 

(|^(0)| 2 ,...,|^(0)| 2 )eS m (28) 
(|Vi(T)| 2 ,...,|^(T)| 2 ) eS fm (29) 

where: 
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• f° depends only on the moduli of controls, 

• there hold (HI'), (H2), (H3). where (HI') is the following condition: 

(HI') the matrix H(t) is skew-Hermitian , measurable as function oft and its elements are either iden- 
tically zero or controls. 

We have the following: 

Theorem 1 (answer yes to Ql: existence of resonant minimizers) Let (ip(.), H(.)), defined in [0, T], 
be a solution to the minimization problem (P'). Then there exists another solution (tp(.), H(.)), to the same 
minimization problem (with the same source, target and initial condition tp(0) = tp(0)) that is in resonance: 

H j ,k(t) = U j , k (t)e i ^ k , where: (30) 
U jtk (.) : [0, T] R, Uj, k = -U kJ (31) 
<Pj, k := arg{^{0)) - arg(^ k {Q)) G [-n,ir] (32) 

Moreover the arguments oftpj{.) do not depend on the time. In formula (32), ifipj(0) = 0, then arg(ipj(0)) is 
intended to be an arbitrary number. 

Theorem 2 (answer to Q2) Let (ip(. ), H(.)), defined in [0, T], be a solution to the minimization problem 
(P'). Then if the cost is strictly increasing, (tp(.), H(.)) is weakly-resonant: 

H 3At)\i jihil =Uj,k,l(t> iVi - k -\ where: (33) 
U 3 , k ,i{.) : Ij, k ,i -> R, U j>k> i = -U k , jt i (34) 
<Pj,k,l ■= arg(ipj(a jtkt i)) - arg(ip k (a,j, k ,i)) G [-tt,tt] (35) 

where Ij, k i =]dj k,U^j kl[ are ^ e intervals defined in Section 1.5.1. In formula (35), if ipj(&j,k,i) — 0, then 
arg((tl)j{a,j t k,i))) is intended to be an arbitrary number. 

In the rest of the paper we call an admissible pair, a couple trajectory-control for an optimal control problem 
satisfying (H1'),(H2),(H3). 

Remark 15 Notice that after elimination of the drift, to be in resonance implies to use controls with constant 
phases. 

An important consequence of Theorem 1 is the following. If an admissible pair is a solution in resonance to 
the minimization problem, then there exists also one for which axg(ipj(0)) = 0, j = 1, ...,n, i.e. corresponding 
to real controls and real initial conditions . In this case equation (26) restricts to reals, so ip(t) € S*" -1 C R n . 
Hence it makes sense to give the following: 

Definition 3 We call (RP') ; the problem (P') in which all the coordinates and controls are real. 
And we have: 

Corollary 1 If (?/>(.), H(.)) is a solution to the minimization problem (RP') (that is with ip(f) £ S n ~ 1 G 
R", H(t) G so{n)) then it is also a solution to the original problem (P') or equivalently, with the transformation 
described in Section 2, of the original problem (P). 

The following Theorem restricts the set of candidates strictly abnormal minimizers, for (RP'): 

Theorem 3 (partial answer to Q3 for (RP')) Let {ip(.),H(.)) be a couple trajectory-control that is a 
minimizer for (RP') and assume that there are no constraints on the controls (H(.)j tk : Domfy) — > RJ. Then 
for every t G Dom(?p), there exists an interval [t\, £2] arbitrarily close to t (possibly non containing t), on which 
("0(-)i H (.)) is not a stricly abnormal extremal. 

For properties of abnormal extremals, we refer the reader to [1, 8, 9, 27, 28]. 
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Remark 16 In Theorem 3 the fact that we assume unconstrained controls is just a technical hypothesis that 
permits to simplify statements and proofs. Notice that for the costs described in Section 1.4 controls can always 
be assumed without constraints. 

As we said above, the non existence of stricly abnormal minimizer is still a conjecture. After the result given 
by Theorem 3 to prove this conjecture one is essentially left with the following: 

Problem: Let (ip(-), H(.)) and (tp(.), H(.)) be two non strictly abnormal minimizers defined respectively 
on [i a ,ib] and [ib,i c ] with t a < tb < t c . Is it true that the concatenation of the two is not strictly abnormal? 



4 Resonance and Weak- resonance 

In this Section we prove Theorems 1 and 2. The key point is to identify the components of the control Hj t k(-) 
that are responsible of the evolution of and of &rg(ipj). The difficulty is that these components are well 
defined only for the times such that 7^ an d i^k{t) ^ (ipj and tpk are the states coupled by H^k{-))- So 

we have to split the problem into the intervals Ij t k,i defined in Section 1.5.1. 

Let (tp(.), H(.)) be a solution to the minimization problem and restrict the problem to the interval I a , 
a = (j, fc, I) £ (1, ...n) 2 x (1, to), where ipj and ip k never vanish. Define in I a : 

6f(t) := argty?(t)), 
p?(t) := |^(*)|, 
/J&(t) :=9f(t)-em- 

From the definition of I a the phases arc well defined. Define the functions u" k (t) and v" k (t) in such a way 
that: 



We have: 



Now in I n 



Hf >k {t) := H jtk (t)L =: («? fc (t) + i^ k {t))e*W\ (36) 



U j,k - U k,j> V j,k~ V k,j fo 7 \ 

\ H r, k \ = Kk+Kk\- ( ' 



Mt) =Y, H Ut)Mt) = E + ivf tk {t))\Mt)\^ (t) - 

k k 

From that, with a simple computation one gets: 

k 
k 

These equations show that, on I a , u^ k and v^ k are responsible respectively of the evolution of \ipj | and arg{ipj). 

From the original minimizer (ip(.), H(.)) we now define a new minimizer by means of two consecutive 
transformations: 

• the first in which we set at zero all the ^"-components of the control and we get a new admissible pair 
connecting Si n and Sfi n having a strictly smaller cost for a strictly increasing / and that is weakly- 
resonant. This proves Theorem 2; 

• the second in which Hj tk is set to zero where ipjipk = and in which all the phases /3" fe are set to 
arg(ipj(0)) — arg(ijj k {^))- This transformation does not increase the cost and produce a control in reso- 
nance. This proves Theorem 1. 
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These transformations are realized in the following. 

Let (V'(-)j be the original minimizer and let u" k and v" k the corresponding components of the control in 
I a , defined in formula (36). Define a new admissible pair (■tp(.), H(.)) by setting in I a : 

Vf tk = [W> 

Notice the important point that v" k = implies that the corresponding (3" k are constants (see equation (39)) 
and they will be considered arbitrary phases. In other words (ip(-), H (.)) is the admissible pair corresponding 
to the control: 

]M ' 1 - a i (t)e i ^ iitel a for some a, 



where j3" k are constant arbitrary phases. The corresponding trajectory in each I a satisfies: 

r u*)=m)<& (4 2 ) 

with and 6 k constant and such that 6" — 8 k = (3" k . Now since Hj t k = Hj.k if t ^ U a I a and the Uj^ 
component is always the same for every a, it follows that |V>j(t)| = |^(i)| for every t G [0, Tj. This means that 
ip(.) and %p(.) connect the same source and target (recall that they are defined by conditions on moduli). By 
construction (■)/>(.), if (.)) is in weak-resonance . 

Proof of Theorem 2 By contradiction assume that: 
• /° is strictly increasing; 



(ip(-), H (.)) is not in weak-resonance that means v" k {t) ^ on some set of positive measure A C [0,T] 
and for some indexes j, k, a. 



In this case from equation (37) we have: 

\Rj,k(t)\ < \H jtk (t)\ for each t G A, 

that means: 

f(H(t))dt< [ f(H(t))dt, 



io Jo 

which contradicts the optimality of (4>(-), H(.)), proving Theorem 2. H 

Proof of Theorem 1 In formula (41), the phases [3" k are arbitrary so for each interval they could be 
chosen to be: 

Pj,k = <Pi,k ■= arg{iPM) ~ arg(M<>)) e [-%, tt] (43) 

where arg(ijjj (0)) is intended to be an arbitrary number if ipj(0) = 0. This proves that if (tp(.), H(.)) is a 
solution to the minimization problem, then there exists another solution (vKOj that correspond to control 

H of the form (30), (31), (32) on \J a I a . It remains to consider the set Badj,k =]0, T[\ U a I a , already defined in 
formula (17), where tpj or ip^ are zero. The following Lemma assures that on Badj^, we can put the control 
at zero, without changing ip and without increasing the cost. 

Lemma 1 Let (ip(. ), H(.)) be an admissible pair and define a new control H(.) by: 

/y , ( Hj^{t) for all t such that ipj(t)ipk(t) ^ 

3 ' ' | otherwise. 

Then (ip(-), H(.)) is also an admissible pair and its cost is not bigger than the cost of (tp(.), H(.)). 
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Proof. Let us first prove that H(.) is a measurable function. Let Nj be the set of times where tpj is not 0. Being 
ipj(-) measurable it follows that the function 1L, defined by TLj(t) = 1 if t G Nj, TLj(t) — if t £ Nj, is measur- 
able. This implies that Hj ik (t) = ^j(t)^k(t)Hj^ k (t) is measurable and H(t) as well. Let us now recall a classical: 

Fact. Let g(.) : M — > M be an absolutely continuous function. Then the set of points x where g(x) = 
and g'(x) ^ has zero measure. 

From this fact, we have a. e. that ^ fe Hj, k (t)ip k (t) = ^2 k Hj tk {t)ij} k {t). Indeed, when ifjk(t) = then Hj^ k (t)ip k (t) - 
Hj,k(t)ijjk(t) = and for a.e. t such that ipj(t) = we have ipj(t) = and so ^2 k Hj^ k (t)ijjk(t) = = 
J2kHj,k{t)iJk{t). Hence H(t)ip(t) = H(t)ip(t) for a. c. t in [0,T]. It follows that if ip(t) = H(t)ip(t) a.e., then 
this relation holds also with H(t) in the place of H(t), that is (ip(-), H(.)) is an admissible pair. From the fact 
that \Hj k \ < \Hj k\, it follows that H(.)) has a cost not bigger than (ip(.), H(.)). The Lemma is proved. □ 

The proof of Theorem 1 is now complete. 



Remark 17 In formula (43), /3" fe has been chosen in such a way to be compatible with the initial condition: 

^(0) = V(0). (44) 

4.1 Trans versality 

In this Section we give a more geometric interpretation to the proofs of Theorems 1 and 2. We are going to show 
that the vector field associated to the control v^ k defined in formula (36) is always tangent to a submanifold 
of S 2n ~ 1 whose points are reached with the same cost. As a consequence the transversality and maximum 
conditions of PMP say that the control v" k can (resp. must) be set to zero for a nondecreasing (resp strictly 
increasing) cost. 

Lemma 2 Given cti,...,a n S [— tt,tt], let us define the map: 

Kota : { (^i,...,^„) ~ (e^Vi,-^^)- 

If ip{.) = (tpi (.),..., tpn{-)) defined on [0, T] is an admissible curve then Rot a (ip(.)) is also an admissible curve 
and has the same cost as ^(.). 

Proof. It is just a matter of computation. Let us denote by H (.) the control matrix associated to the admissible 
curve ip(-). The fact that tp(.) is an admissible curve, associated to the controls H(.), writes ip k = J^j Hj^ji 
Vfc e {1, . . . , n}. By multiplying by e iak , we find tpk^ iak = J2j Hj,k^ ak ~ aj ^jG iaj ■ Hence the curve Rot a (ip) 

is an admissible curve corresponding to controls Hj. k = Hj^ k c i ( ak ~ aj \ And since the cost function does not 
depend on the phase of the controls, the two curves have the same cost. 

This means that Rot a is an "isometry" for the cost defined by / , for any (ai, . . . , a n ). A trivial consequence 
of this lemma is: 

Corollary 2 For any two points ip 1 and ip 2 in S" 2 ™ -1 , all the points of the set: 

T r = {(^je' a \. . . ,^ 2 e' Q ") | (ai, . . . ,a n ) G M"} 
are reached with the same cost from the set: 

V ={(Vi^° 1 ,...,^") I («!,..., «n) 6 1"}. 
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As a consequence if ip(.) : [0,T] — > S 2n 1 is a minimizing trajectory between the two sets 7y,i 7"^2, then the 
transversality condition of PMP (see v)) is: 

<P(t),TT m >=0. (45) 

Now let F" k (ip) and G" k (ip) defined in I" k be the two vector fields associated with the controls u" k and v^ k 
defined in formula (36): 

k 

One can easily check that on I a : 

G? fe W) = i(e^^4 + e-^i^)- 

And with a simple computation one can see that the vector Gj k (ip(t)) is tangent to the set T^ t y. 

Gl k m))eTT m , ViG/ a . (46) 

From (45) and (46), one get that < P(t), G? k (ip{t)) >= 0. 

close to any time t of the domain of a given minimizer, there exist an interval of time where the minimizcr 
is not strictly abnormal 

there are not strict abnormal extremals. Then the maximality condition iii) of PMP implies (here / = 
ul k F« k &)+vl k Gl k m-. 

Y, u tk{t) < P(t),F% k Mt)) > +Pof(\uU t )+ iv f,k(t)\) = .max £u i)fc < P(t),F£ k (tP(t)) > +Pof(\u j , k +iv j , 

j,k j,k 

Now since po < (there are no strictly abnormal extremal, as it will be proved in Section 5), we get that (we 
stress the fact that v^ k is not useful to reach the final target since the source and the target are defined by 
conditions on the moduli and G" k is responsible only of the evolution of the phases, see eq (39)): 

• for a non decreasing cost this condition is realized if v" k = a.e. (plus conditions on u" k ); 

• for a strictly increasing cost this condition can be realized only if u" fc = a.e. 

4.2 Eigenstates 

Let us recall that an eigenstate is a state for which one of the coordinates has norm 1 and the others are 0. 

In this Section we show that it is possible to join every couple of eigenstates by a minimizing trajectory 
that is in resonance (question Ql'). More in general we prove that the answer to question Ql' is yes for any 
couple of initial and final states ip 1 = (ipl, ip^) and ip 2 = (V'ij sucn that ^Pj^Pj = f° r an y Jj that in 

particular is true for eigenstates. 

Let tp(.) be a minimizing trajectory in resonance between the points (\ip\ |, IV^I) an d (l^i |> •••) IVnD- Define 
9j and 9 2 as the arguments of ipj and ip 2 , putting them to when the corresponding coordinate is 0. We define 
otj to be equal to 0j if 0j ^ 0, or to 9 2 if ip? ^ and if both tpj and tp 2 are 0. Finally let a := {a\, a n ). 
The curve 

</>(•) := Rot a $(.)), 

is a resonant minimizer for the problem with initial and final condition ip 1 and ip 2 . 

4.3 An Example of a Non- Weakly-Resonant Minimizer 

In this Section, in the case of a non strictly increasing cost, we show an example of an optimal couple trajectory- 
control (ip(-), H[.)), joining a source and a target defined by conditions on the moduli, that is neither resonant 
nor weakly-resonant. Consider a time minimization problem for a 4-level system in the form (2) (i.e. with 
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controls on the lower and upper diagonal), for the isotropic case. We have f°{H) = max-fjifi^l, 1^2,3 1, |^3,4|} 5 
that is not a strictly increasing function but just a not decreasing one. 

It is easy to see that the following curve: t t— > (cos(t), sin(t), 0, 0) for t G [0, |] is a minimizer between 
the eigenstates Vi| = 1 ancl |V'2| = 1. It can be obtained by different control functions. For instance by 
t i ► (-1,0,0), or by 1 1-> (-1,0, f/ 3 ) where 



The second one is clearly not weakly-resonant. 

The following is still an open question: Under the hypothesis (HI'), (H2) and (H3), does there exist an 
example of non-decreasing cost function such that there exists a trajectory H(.)), solution of (P'), being 

non- weakly-resonant, such that there is no weakly-resonant control H with (if)(.), H(.)) solution of (P') ? 

5 Strictly Abnormal Minimizers for the real problem 

In this Section we prove Theorem 3 (i.e. we prove that close to any time t of the domain of a given minimizer, 
there exist an interval [<i,i2] where the minimizer is not strictly abnormal). We are able to prove the result 
only in such a interval, since on that interval we can make a suitable partition of indexes, see Lemma 3 (how to 
extend this result to the whole domain of the minimizer is still an open question, cf. end of Section 3). Then 
the difficulty of the proof is coming from the fact that the dynamics has singularities each time a coordinate 
is zero. For instance, in the most important example (see formula (2)), when ipi = V>2 = then the control 
Vi t 2 has no effect on the dynamics, i.e. the corresponding vector field is zero when ipi = ip2 = 0. Hence, in 
order to prove the theorem, we show (subsections 5.1 and 5.2) that every minimizer tp of the problem (RP') 
is also a solution of an auxiliary optimal problem (RP"), living on a submanifold J of S" 1-1 , which has no 
singularities. Then in subsection 5.3 we prove that this new problem has no abnormal extremals (which proves 
that the minimizer ip is a normal extremal as solution of (RP")). Finally, we come back to the original problem 
(RP') proving that the minimizer has a normal lift in T*S n ~ 1 (subsection 5.4). 

In the following, let H(.)) be a couple trajectory-control that is a minimizer for (RP'). 

5.1 Permutation of indexes 

In this subsection, we construct a submanifold J on which we will restrict the optimal problem (RP'). 

Lemma 3 In every neighborhood of every time t of the domain of the minimizer, there exists a sub-interval 
of time, denoted by [ti,t2] (possibly non containing t), such that there exists a partition I {J J of {l,...,n} 
satisfying: 

• if j G / then ipj{t) = Vt G [ti,t 2 ], 

• if j G J then ipj(t) Vt G [h,t 2 ]. 

This is just a consequence of the continuity of The proof is left to the reader. In the following we always 
restrict to the interval [ti,^]- 

Definition 4 We say that two indexes j and k of J are connected ( denoted by j ~ k) if j = k or there exists 
a sequence ji, j s of indexes of J such that j = j\, k = j s and Vr < s the coefficient Hj T j r+1 is a control. 

Remark 18 This definition of connectedness is exactly the same as in subsection 1.2 but for the sub-graph 
defined by J. Until the end of this Section, "connected" refers to this last definition. 

We immediately get: 



u 3 (t) 
u 3 (t) 

W) 
u 3 (t) 



I 



1 



1 



forte [0,f], 
forte [f,f], 
fortG[f,f], 
fortG[f,f]. 
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Lemma 4 In J, the relation "<~ " is an equivalence relation. 

Definition 5 We denote K\,...,K r the equivalence classes defined by and mi, ...,m r their cardinalities. 
We also denote Mq = and = X)fc</ TO fc- 

Now, let us make a permutation on the indexes that orders the sets K\, K r , I in such a way that: 

W<r K^fM^ + l;...;^}, 

and 

7 = {M r + l;...;n}. 

Lemma 5 After the permutation, we have: 

a) /or aZZ j > M r + 1 (j S I), V>j = on [*i,i 2 ]. 

b) /or a// i < r, the map is constant on [ii,i 2 ]. 

j£K e j=M e _ 1 + l 

Proof. The point a) is just a consequence of the definition of I. In order to prove the point b), let us consider 
£ < r. For any j in Ki, we have: 

k<n 

But Hj t k = if k £ Ki U I and ^ = if k € I. Hence for t in [ti, i 2 ]: 

few = X H^ k (t)Mt). 

keK t 

Now, the matrix H l (t), whose coefficients are the Hj^(t) with j, fc £ iQ, belongs to so{mi). Hence the vec- 
tor ijj , whose coefficients are the i/'jW with j E iff, satisfies ^ = H e ip e and then has constant norm, i.e. 
Sje_fs" f hAj'W | 2 i s constant. H 

From the previous Lemma 5, it follows that for t G [ii, i 2 ], belongs to the set: 

J = r 1 " 1 ^) X ... X S^-^Cr) X HWjih)}, 

where Cfc = yjT,jeK t l^(*i)| 2 - 

Remark 19 In the following, we are going to restrict (RP') to J . Notice that the dimension of J is X^<r( m ^ — 
1) < n — 1, while the original problem (RP') lives on S" -1 that has dimension n — 1. 

5.2 Restriction of the Control System to 

In this Section, we show that ip is also a solution of an auxiliary optimal problem (RP") in J . 

Definition 6 We call (RP") the optimal problem (RP') in which we restrict the dynamics by adding a new 
condition on the matrix of controls: Hjj.(t) is set to if j or k is in I . 

Remark 20 Notice that J is preserved by the dynamics of (RP"). 

Lemma 6 The curve is a minimizes for (RP"). 

Proof. As proved in Lemma 1, ip does not change if we set to zero the controls flj - ,fc such that one index 
belongs to I. Hence ip is an admissible curve for (RP"). 

Now, since the cost function is the same for both (RP') and (RP"), it follows that ip is a minimizer for 
(RP"). 
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5.3 The Minimizer ^(.) is not an Abnormal Extremal for (RP") 

Let us denote Fj^{ip) the vector field associated with the control Hj k- 4> = Y^j k Hj,kFj,k(' l P)- We denote by 
Ag the distribution generated by the Fj^ip) with j, k G Kg and A = @iAg_. We have that A(q) C T q J for 

gej. 

Using Fj^(if>), the dynamics of (RP") reads: 

£ j,k<EK e 

Lemma 7 The distribution A is £/ie whole tangent space of J : A = T J . 
Proof. Let us first prove that, for every £, A £ = TS™*- 1 ^). 

First step: it exists a set of couples of indexes of Kg, with cardinality mg — 1, that keeps -K^ connected in 
the following sense: for any couple of indexes j, k in -fQ, there is a sequence ji, j s in i£t such that j\ = j, 
3s = k, {jp,j P +i} 6 i ; and #j Wp +1 is a control. 

The proof can be done by induction on the cardinality of Kg. If it has cardinality one. Lg is empty. If it has 
cardinality two, the proof is trivial. If the cardinality of Kg is bigger than two, we can choose an index k G Kg 
such that Kg := Kg — {k} is still connected (this is a standard fact of graph theory). Now, Kg is connected 
and has cardinality mg — 1. Hence, by induction hypothesis, there exists an Lg that keeps Kg connected and 
has cardinality mg — 2. Now, k is connected to an index j of Kg, hence the set Lg = Lg U {{j, k}} keeps Kg 
connected and has the required cardinality. 

Second step: the family of Fj^itp) ({j> ^} £ Lt) is linearly independent (Lg was constructed for this purpose). 

The proof can be done by induction : if mg = 1, Lg is empty. If mg > 1, then one can choose a index k 
appearing only once in Lg. Let call j the index such that {j, k} G Lg. The vector Fj t k(ifi) is linearly independent 
of the rest of the family : it is the only one vector that has a non zero /c-coordinate. 

Now Lg — {{j, k}} has cardinality mg — 2 and keeps connected Kg — {k} of cardinality mg — 1. Hence we 
can apply induction. 

Last step: From the previous step, it follows that Ag has the same cardinality mg — 1 as the sphere S me ~ 1 (Cg), 
hence we get Ag = TS" 1 ^ 1 ^). 

Thus we have A = ®gAg = (SgTS" 1 '^ 1 ^) = TJ. ■ 

Since the distribution A is the whole tangent space of J , it is a standard fact that there are no abnormal 
extremals for the problem (RP"). 

5.4 End of the proof 

Now, we are going to prove that the curve ip(.) admits a normal lift on S*™ -1 i.e. it is not a strictly abnormal 
extremal for the problem (RP'). 

Let us denote a local coordinate system on S me ~ 1 (Cg), and O = (Oi, r ). Then: 

(6i, r , Cl, C r , 1pM r+1 , V'n) 

is a local coordinate system on K™. We denote by Pq^Pci and P^,. the dual coordinates in T*M. n , 

l<r £<r i£l 

and 

Pj = Y, p B t dQi 

i<r 

its restriction to J . 
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Remark 21 Notice that the vector fields Fj^ with j, k in Kg., depend only on the 0f coordinates and that they 
annihilate dC s (s < r) and dipi (i S I)- This fact will be essential to conclude the proof. 

The PMP-Hamiltonian associated with (RP") is: 

Mj(Q,H,Pj,p ) = Pj( H hk F hk {Q))+ Po f(H), 

j,ke.J 

where /° is the cost function. Since ip is a normal extremal for (RP"), there is a lift (t/j, H, Pj,po), with 
Po ^ 0, that satisfies: 



Pe e = -^(e,H,Pj,p ) 



and the maximahty condition: 



Mj(ip(t), H(t), Pj(t),p ) = max{M^(V(i), H, Pj(t),p )}, 

H 



where the maximum is taken over the set of controls. 
Now, the PMP-Hamiltonian associated with (RP') is: 

m(tP,H,P,p ) =P(Y^ H 3,kFj,k) +Pof°(H). 

j,k<n 

To conclude the proof of Theorem 3 we need the following: 

Claim Let us denote: P = Pj + Yli<r dC( + Yliei dipi. Then [ij),H, P,po) is a normal lift of ip sat- 
isfying the PMP for (RP'). 

Proof. Let us first remark that P satisfies Pc e = P^ t = for every £ < r and i € I. 

In order to prove the claim, we have to prove that the lift satisfies the Hamiltonian equations: 

$ = ™&n,P,Po) (47) 

P = ~(if,H,P,Po) (48) 

and the maximality condition: 

m$(t),H(t),P(t),p ) =irmx{m( : >P(t),H,P(t),po)}. (49) 

H 

The equation (47) is always trivially satisfied (it is the dynamics). 
Let us prove that (48) holds. Firstly, we have that: 

K = -^(Q,H,Pj,p ) = -^$ t H,P,po). 

Secondly, 

j,k<n 



but Hj k = if i or k is in / and, because of Remark 21, i^r(Fj,k) = when j and k arc in J. Hence 

oc, 



Pc e = = —§ M -(^P, H, P,po)- The proof is the same for P^ li for i & I. Hence (48) is satisfied 
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Let us now prove (49). We have that, if j, or k, is in I then Fj^{ip) is in span{g^- ; i E I}. Indeed, 
Fj,k$) = $k-$r ~ ^af^. and tlicn if .? e 7 - then = and FjA^) = ^fcafj- Hence, because = 0, we 
have that P(Fj t k) = if j or fc is in / and we have that for any control H: 

j,k<n 

= P(J2 H^F^))+p f (H) 
= Py( J2 H hk F. hk ^))+p Q f(H) 
= Mj(ij,H,Pj,po) 

Now the fact that (V 1 , H, Pj,Po) satisfies the maximality condition for (RP") allows to conclude that (ifi, H, P,po) 
satisfies the maximality condition for (RP'). H 
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